Ferroelectric Phase Transitions in Ultra-thin Films of BaTiOs 
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We present molecular dynamics simulations of a realistic model of an ultrathin film of BaTiOa 
sandwiched between short-circuited electrodes to determine and understand effects of film thickness, 
epitaxial strain and the nature of electrodes on its ferroelectric phase transitions as a function of 
temperature. We determine a full epitaxial strain-temperature phase diagram in the presence of 
perfect electrodes. Even with the vanishing depolarization field, we find that ferroelectric phase 
transitions to states with in-plane and out-of-plane components of polarization exhibit dependence 
on thickness; it arises from the interactions of local dipoles with their electrostatic images in the 
presence of electrodes. Secondly, in the presence of relatively bad metal electrodes which only 
partly compensate the surface charges and depolarization field, a qualitatively different phase with 
stripe-like domains is stabilized at low temperature. 
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Ferroelectric (FE) materials are vital to the technolo- 
gies based on micro and nano electromechanical systems 
and certain random access memories. The fact that fer- 
roelectricity in finite objects is very sensitive to mechan- 
ical and electrical boundary conditions is of great rele- 
vance to nano-scale devices based on ferroelectric mate- 
rials [1, 2j. These boundary conditions in the case of a 
thin film are determined by the properties of its interfaces 
with substrate on which it is grown and the electrodes. 
While understanding the suppression of ferroelectricity 
in perovskite thin films has been a fundamental issue [3| , 
a clever choice of an electrode and a substrate can be 
used effectively in designing nano-scale ferroelectric film 
based structures with desired properties. For example, a 
high-Tc lead-free ferroelectric was developed by growing 
BaTiOa films on an appropriate substrate Q. 

FE materials are spontaneously polarized below a cer- 
tain temperature and the direction of this spontaneous 
polarization can be switched by an external electric field 
[5|. When a FE film is polarized in a direction perpen- 
dicular to its plane, the bound charges at the surface of 
the film give rise to depolarization field that suppresses 
the spontaneous polarization. Free carriers in the in- 
terfacing electrodes partially compensate these surface 
charges, and aid in maintaining the ferroelectricity down 
to 6 unit cell thicknesses [6,]- While the suppression of 
the out-of-plane polarization has been investigated fairly 
well [1, 0; Q and can be understood with relative ease, 
effects of electrodes on the in-plane polarization and asso- 
ciated ferroelectricity have not been explored yet; their 
origin is expected to be fundamentally different. Since 
the switching of polarization necessary for memory ap- 
plications occurs by rotation Q through the states with 
in-plane polarization, understanding properties of these 



states is very important to applications of FE materials. 

BaTiOa, a standard example of FE perovskite ex- 
hibits a sequence of phase transitions from cubic para- 
electric to tetragonal (C<->T) to orthorhombic (T<- 
>0) to rhombohedral (0<->R) ferroelectric phases 
characterized by polarization ordering along (001), (110) 
and (111) directions respectively. It is an ideal case 
for exploration of effects of electrodes, epitaxial strain 
and film thickness on in-plane ferroelectricity. Phe- 
nomenological Landau theory has been used to map 
a temperature-misfit strain phase diagram of BaTiOs 
and PbTiOs films of thickness > 50 nm grown on cu- 
bic substrate, with short circuit boundary conditions 
Q. Dieguez et al., mapped the equilibrium BaTiOa 
structure at zero temperature as a function of epitax- 
ial strain and simulated corresponding ferroelectric tran- 
sitions in the bulk using the effective-Hamiltonian ap- 
proach [l^, |Tl|, [l^. Information about the strain- 
temperature phase diagram is very useful in the design 
of ferroelectric nano-structures. 

In this letter, we use a first-principles model Hamil- 
tonian of BaTiOs of Ref. [13 along with a simple elec- 
trostatic model of electrodes [ij] in classical molecular 
dynamics simulations and determine (a) thickness depen- 
dence of ferroelectric transitions in BaTiOs films sand- 
wiched between limiting cases of electrodes (perfect vs. 
imperfect), (b) a complete epitaxial strain-temperature 
phase diagram in the case of perfect electrodes and (c) 
the nature of low-temperature phase and ferroelectric- 
ity when the electrodes are not so perfect. We point out 
that interactions with electrostatic images of the in-plane 
dipoles interacting are responsible for size dependent fer- 
roelectricity, even in the absence of depolarization field. 

We used an effective lattice dynamical (phonon) 
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FIG. 1: Schematic representation of a FE thin film of thick- 
ness I unit cells (here Z=2) sandwiched between (a) perfect and 
(b) imperfect electrodes. Horizontal thick lines marked as "E" 
represent the electrostatic mirrors used to model electrodes. 
They are distance | away from the FE film surface (d=0 in 
(a), d=\ in (b)). Each arrow represents a local dipole within 
a unit cell (a=3.94 A ) of the FE crystal. Thick dashed lines 
indicate boundaries of the periodic box of the system used in 
simulations. 



Hamiltonian of BaTiOs expressed in terms of local- 
ized atomic displacement basis spanned by two classi- 
cal 3-dimensional vector degrees of freedom (DoF) per 
unit cell: one representing the unstable polar optical 
branches of the cubic structure and another represent- 
ing the branches of acoustic phonons. This amounts to 
approximately integrating out the high energy phonons, 
reducing DoFs from 15 to 6 per unit cell. The average 
value of the former yields polarization, where as the lat- 
ter describes local strains in the crystal. In addition, 
six components of the homogeneous strain tensor are in- 
cluded in the DoFs. This Hamiltonian includes harmonic 
and anharmonic interactions of phonons, coupling be- 
tween polar and strain DoFs, and elastic contribution to 
energy. We use the parameters in effective Hamiltonian 
as determined earlier from first-principles local density 
functional theory calculations [lO, [Hi . Due to underes- 
timation of lattice constant in this scheme, a negative 
pressure of -5 GPa is used in all simulations for compar- 
ison with experiment. An effective mass is attached to 
each polar DoF, based on the eigenvector of the unstable 
mode at F— point and ionic masses. 

To model a FE thin film sandwiched between two elec- 
trodes, we have treated electrodes as perfect electrostatic 
mirrors located at a distance | away from the surface of 
the film (see Fig. 1). This "gap" can also be thought of as 
a dead layer separating electrodes from the FE film. The 
length-scale d determines how effective the electrodes are 
in compensating the depolarization field p^ : 
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where Pz is the out-of-plane polarization in the film, I 
being the thickness of film (both I and d are in the units 
of number of unit cells). Here, we use two limiting values 
of d: d — corresponding to "perfect" electrodes and 
vanishing depolarization field, and d = I corresponding 
to "bad" electrodes on the BaTiOa thin films. It is clear 
that an electrostatic image of an electric dipole in the 
FE film is a dipole with same magnitude inside the elec- 
trode. However, the direction of the image dipole is the 
same (opposite) as that of the dipole in the film, if it is in 
z (x or y) direction. In this work, we consider identical 
electrodes on the two sides of the FE film, which allows 
simulations with periodic boundary conditions (due to 
infinite number of images, see Fig.l) with a periodicity 
of 2(l+d), i.e., treatment of electrodes amounts to simu- 
lating periodic system with twice the size of the film. We 
simulated FE films of two types: (a) bulk-like (system F) 
which have no epitaxial constraint and can be strained 
in the plane of the films, and (b) epitaxial films (sys- 
tem EF), whose in-plane strain is fixed by the choice of 
the substrate (electrode). We keep the in-plane homoge- 
neous strain DoF fixed during the simulation of the latter. 
As a check on internal consistency and for comparison, 
we also simulated ferroelectricity in bulk BaTiOa (sys- 
tem B) with the same size(s) using periodic {not mirror) 
boundary conditions (note that there are no electrodes 
in system B, unlike in systems F and EF). 

Mixed-space molecular dynamics, used earlier in large- 
scale simulations of relaxors [3, [lB| treating infinite- 
range dipolar interactions in reciprocal space, is used 
here to determine finite temperature properties of the ef- 
fective Hamiltonian. The Nose-Poincare fi3\ thermostat 
with simplectic properties implemented in the present 
simulations allows us to use a relatively large timestep 
of 2 femtoseconds. At each temperature, the system 
was thermalized with 50,000 timesteps and averaging was 
performed using configurations of subsequent 150,000 
timesteps (amounting to a simulation of 0.3 nano sec). 
Since ferroelectric phase transitions are known to be first- 
order and exhibit hysteresis, we increased (decreased) 
temperature in small steps of 10 to 20 K in the heat- 
ing (cooling) runs of simulations, with smaller tempera- 
ture steps near the transition. The size of systems in our 
simulations is 16 x 16 x /, where I = 2, 11 unit cells. 

We first report our results for FE films sandwiched be- 
tween perfect electrodes: d — 0. In the limit of large 
thickness, ferroelectric properties of thin films with no 
epitaxial constraint (system F) are expected to converge 
to the bulk behavior. We find for both the systems F 
and B that the transition temperatures do depend on 
the thickness for polarization ordering in all the three 
directions of polarization (Fig. 2a). The dependence of 
Tc of the system B arises from the statistical mechani- 
cal effects of the finite size [l3|, where as the difference 
between the Tc's of systems B and F reveal the nano- 
size effects of the film thickness. It is evident that the 
transitions to FE states with in-plane polarization occur- 
ring at higher temperatures in the films converge to bulk 
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FIG. 2: Dependence of transition temperatures on the thick- 
ness of films sandwiched by perfect electrodes (a) with no and 
(b) with epitaxial constraint (B: Bulk, given for comparison, 
F: Film without epitaxial strain, EF: Films with epitaxial 
strain): a = lattice constant of BaTiOa. 



behavior when the film thickness is greater than 10 unit 
cell layers, where as the transition to the state with out- 
of-plane polarization converges more gradually to that in 
the bulk. This dependence of Tdz) on the film thick- 
ness even in the case of vanishing depolarization field is 
surprising. 

A related finding from our simulations is that even the 
high temperature paraelectric phase P = of films F 
exhibits tetragonality through uniaxial symmetry of its 
dielectric constant. The origin of both the findings can be 
traced to the fact that electrostatic images of the in-plane 
local dipoles in the presence of electrodes are inverted and 
break the horizontal reflection symmetry ah ■ Since these 
images interact with local dipoles in the film polarized 
in any general direction, in a way distinct from that in 
the bulk, the presence of electrodes gives rise to thickness 
dependence of ferroelectric transitions in films even in the 
absence of depolarization field. 

We find a rather different thickness dependence of the 
ferroelectric transition behavior in epitaxial BaTiOs films 
sandwiched between perfect electrodes (see Fig. 2b). In 
this case, the in-plane strain is clamped which suppresses 
the fluctuations associated with in-plane dipole moments 
(FE transitions are fluctuation driven first-order phase 
transitions [l^l) and hence the transition temperature for 
in-plane ordering is reduced, which remains almost con- 
stant at 70 K (±10 K) as a function of film thicknesses. 
Unlike bulk, there is no orthorhombic phase stable at in- 
termediate temperatures. Epitaxial constraint naturally 
favors uniaxiality, resulting in a noticeable enhancement 
of the temperature of the transition (well above the bulk 
transition) to the state with out-of-plane polarization, 
consistent with experimental observations Our sim- 
ulations suggest a critical thickness of 3 unit cell layers, 
below which we find no ordering of Pz at low tempera- 



tures. 

The value of epitaxial strain has a spectacular effect on 
ferroelectric phase transitions in the films (see Fig. 3 for 
results in the presence of perfect electrodes d = 0). At 
compressive epitaxial strains {exx < —0.005), BaTiOs 
films exhibit a single transition to tetragonal phase, 
with increasing Tc with the magnitude of strain. For 
—0.005 < exx < 0.03, they exhibit two transitions from 
paraelectric to tetragonal FE and then distorted rhom- 
bohedral FE phases (i? ) as the temperature is lowered. 
There is a crossover point at exx = 0.012, below (above) 
which the low temperature phase has Pz > Px ^ Py 
{Pz < Px = Py); at this point, there is only one second 
order transition from cubic to rhombohedral FE phase, 
which does not seem to depend on the film thickness. 
For exx > 0.03, there is a single phase transition to or- 
thorhombic FE phase whose Tc increases with the mag- 
nitude of strain. Through comparison with the work 
of Dieguez et al [9i], we show (a) the presence of per- 
fect electrodes results in a shift in the temperature-strain 
phase diagram to the right along strain axis by roughly 
^^xx ~ 0.012, and (b) and there is an overall reduction 
in the transition temperatures by about 70 K partly due 
to the dependence on film thickness. General features 
of our phase diagram are closer to that determined by 
Dieguez et al [9] than the one of Pertsev et al We 
note that the low-temperature aspects of this phase di- 
agram are consistent with properties of BaTiOa/SrTiOs 
superlattice as a function of epitaxial strain [181] . 
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FIG. 3: Temperature-epitaxial strain phase diagram of a film 
(of 4 unit cell thickness) sandwiched between perfect elec- 
trodes (d = 0). Region I is tetragonal, II and III are the R 
phase (see text and [Ref.8]), and IV is orthorhombic. Regions 
II (P^ (= Py) < Pz) and III (P^ (= Py) > Pz) are separated 
by a rhombohedral phase R. 

Ferroelectricity in the ultra-thin films of BaTiOa sand- 
wiched between "bad" electrodes (c? = 1) is qualitatively 
different from its bulk behavior. In this case, the depolar- 
ization field is quite sizeable and effective in suppressing 
Pz- The dielectric constant (see Fig. 4, top panel) e^^ 
exhibits a divergence near T = K, but no signature of a 
phase transition to a ferroelectric state with polarization 
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FIG. 4: Top: Dielectric constant as a function of tempera- 
ture for (left) films (F) and (right) epitaxial films EF {1=4, 
d=l), sandwiched between "bad" electrodes. Bottom: A hor- 
izontal slice on the left {+z and —z polarized sites are sym- 
bolized by empty and filled boxes respectively) and a vertical 
cross-section on the right for 16 x 16 x 4 system for the same 
configuration at 110 K showing domain formation. 



Pz. The tetragonality of the high temperature paraelec- 
tric phase is also apparent in the uniaxial symmetry of 
dielectric constant. The transition to FE phase ordered 
with in-plane polarization occurs at around 70±10 K for 
epitaxial thin films. For the films without epitaxial con- 
straint, this transition occurs at a higher temperature 
210 K). To understand the nature of low temperature 
phase with very high zz dielectric response, we exam- 
ine the snapshots of dipolar configurations of these films. 



They reveal (Fig. 4, bottom panel) formation of vertical 
((010) planes) stripe- like domains with polarization along 
z— axis which form arbitrary patterns in the (001) planes. 
Our results are very similar to those observed experimen- 
tally by Fong et al [Sj. We find that these domains are 
stable throughout the duration of our simulation (about 
0.1 nano sec) and observed at temperatures as high as 
300 K. It is clear that these stripe-like domains are stabi- 
lized through minimization of the energy cost associated 
with depolarization field. Formation of similar domain 
structure was also predicted by Levanyuk et al p^ . 

Origin of most of our observations can be understood 
through Fourier analysis of the dipolar configurations 
noting the mirror boundary conditions and vanishing of 
dipoles in the gap region. For example, Pz can be only 
an even function, where as P^ is always an odd function 
of z: hence q = component of Px has to be zero (for the 
doubly periodic cell). This leads to fluctuations of the in- 
plane electric dipoles whose length-scale is determined by 
the film thickness. Since these dipoles couple with both 
the strain and out-of-plane dipoles, ferroelectricity even 
in the absence of depolarization field exhibits thickness 
dependence. 

In summary, we have shown that the ferroelectricity in 
thin films sandwiched between electrodes exhibits nano- 
thickness dependence even if the depolarization field van- 
ishes. Secondly, the presence of electrodes and their na- 
ture have wide-ranging consequences to ferroelectricity 
in ultra-thin films from simple size dependence to stabi- 
lization of phases with stripe-like domains. The origin of 
these phenomena lies in interaction between local dipoles 
and their electrostatic images in electrodes. 
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